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Abstract. 

We study phonon-assisted dephasing in optically excited semiconductor quantum 
dots within the frameworks of the independent Boson model and optimal control. 
Using a realistic description for the quantum dot states and the phonon coupling, we 
demonstrate that such dephasing has a drastic impact on the coherent optical response. 
We employ optimal control theory to search for control strategies that allow to fight 
decoherence, and show that appropriate tailoring of laser pulses allows a complete 
control of the optical excitation despite the phonon dephasing. 
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1. Introduction 



from lasers [3| over single and entan gled photon sources 



information devices lOj; ll2j; |13|; |l4j; [15 




In semiconductor quantum dots carriers become confined in all three spatial directions, 
resulting in a discrete, atomic-like density of states [H; Q|. For this reason quantum dots 
are often referred to as artificial atoms. Possible applications of quantum dots range 

S 2lB S] to quantum- 
|; |21|. However, while 
for atoms environment couplings can be strongly suppressed by working at ultrahigh 
vacuum and ultralow temperature, for artificial atoms things are more cumbersome 
because they are intimately incorporated in the surrounding solid-state environment 
and suffer from various decoherence channels. This happens even for optical excitation 
of electron-hole pairs in the states of lowest energy, e.g., exciton or biexction groundstate, 
causing the deformation of the surrounding lattice (i.e., formation of a polaron state) 
whereas relaxation is completely inhibited because of the atomic-like carrier density 
of states. In coherent optical spectroscopy 22j; [23|, which is sensitive to the optically 



induced coherence, this partial transfer of quantum coherence from the electron-hole 



state to the lattice degrees of freedom, i.e., phonons, results in dephasing [24]; [25j; [26 



The coupled dot-phonon system is conveniently described within the independent 



Boson model 27|; |28( and possible generalizations 29|; 130). The independent Boson 



32; 33; 34; 35 , whereas 



model is exactly solvable for laser pulses delta-like in time 
approximate description schemes have to be employed for laser pulses of finite duration. 
This was first accomplished within a density- matrix approach by Forstner et al. 36c l37l|. 
who reported a surprisingly large impact of phonon-assisted dephasing on the coherent 
optical response. Apparently, this constitutes a serious drawback for quantum-coherence 
and quantum-information applications in quan tum dots, which have recently received 



considerable interest lid 



l;Q 



18; 19; 21; 38; 39; 40 



A number of quantum control techni ques are known, such as quantum bang-bang 



42 



or spin-echo pulses, 



43 



that allow to 



control, [41( decoherence-free subspaces, 
fight decoherence. However, it is not the system-environment interaction itself that 
leads to decoherence, but the imprint of the quantum state into the environmental 
degrees of freedom: the environment measures the quantum system. Optimal control 



theory 44]; |45|; |46| allows to design control strategies where quantum systems can be 



controlled even in presence of such environment couplings without suffering significant 



decoherence |47|; |48|], e.g. by means of laser of voltage pulse shaping. This surprising 
finding can be attributed to the fact that in the process of decoherence it takes some time 
for the system to become entangled with its environment. If during this entanglement 
buildup the system is acted upon by an appropriately designed control, it becomes 
possible to channel back quantum coherence from the environment to the system. 

It is the purpose of this paper to review the quantum control framework developed 
in refs. [43; |48j and to apply it to optical control of realistic quantum dots. Contrary 
to the simple quantum dot model previously used 43], the strong (bi)exciton-phonon 
coupling will force us to consider not only single- but also multiple-phonon excitations. 
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We have organized our paper as follows. In sec. [2] we briefly review some basics of 
electron-phonon couplings in solids and, more specifically, semiconductor quantum dots, 
and show how to analytically solve the independent Boson model. SectionH]is devoted to 
a description of our approximate wavefunction scheme. We present results for coherent 
optical spectroscopy and compare our results with experiment. Finally, we discuss 
quantum control and optimal quantum control in sec. and summarize in sec. |6j 



2. Theoretical model 

2.1. Electron-phonon interaction 

Electrons in solids experience interactions with the lattice degrees of freedom. Let 

V(r) = Y,V ci (r-R J ) (2.1) 

3 

be the interaction between an electron at position r with the lattice of ions located at 
positions Rj 28]. The potential ^ - V ci (r — Rj) for the periodic lattice, with all ions 



at their equilibrium positions, then gives the bandstructure and the Bloch states of the 
solid. For small displacements of the ions from Rj one can expand the potential (12.11) 
in powers of the ionic displacements Qj, to arrive at the electron-phonon coupling 

V(r) = Y,QrVV ei (r-R j ). (2.2) 

3 

This interaction is to be written in terms of operators. For the phonons we adopt the 
usual displacement operator 

Ql = i J2 (2^)"* (a qjX + al q X ) , (2.3) 

with p the mass density, u q \ the frequency of phonon modes with wavevector q and 
polarization A (longitudinal or transversal), £ q \ the corresponding polarization vectors, 
and a qX the bosonic creation operator. We set h~ = 1 throughout. For the electrons 
and holes of the semiconductor material, we expand the charge fluctuations in terms of 
the fermionic field operators through c k+q c k and d k+q d k , respectively. In the following 
we shall only consider deformation potential coupling to longitudinal acoustic phonons, 



which is known to be the major dephasing channel for carriers in quantum dots 31 



Interactions with longitudinal optical phonons are known to have a profound impact on 



the relaxation dynamics in quantum dots 49|; |50|; l51| . but to be of minor importance 



for the dephasing dynamics of our present concern owing to the flat optical-phonon 
dispersion. The phonon interaction then has the form 

H ep = (2p^)^ \q\ (a, + at,) (d £ c{ +q c k - D h d{ +q d k ) , (2.4) 

where D e and Dh are the deformation potentials for electrons and holes, respectively. 
We have suppressed for conceptual clarity the subscript A in the phonon frequencies 
and operators. Eq. (12.41) has the usual form that either a phonon with wavevector q is 
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Figure 1. Schematic sketch of the excitonic level scheme in a single quantum dot. 
Through proper choice of the light polarization, one can promote the dot from the 
groundstate to either the spin- up or spin-down exciton state X. Within a two- 
photon process one can also directly excite the biexciton state XX of lowest energy, 
whose energy is reduced by the biexciton binding A with respect to twice the exciton 
energy [53 1. 



destroyed and an electron becomes scattered from k to k + qr, or the electron scatters 
from k to k + q upon emitting a phonon with wavevector —q. The same conclusions 
apply to holes. 

2.2. Quantum dot states 

Quantum dots are small islands of lower-bandgap material embedded in a surrounding 



matrix of higher-bandgap material 52|; |53|- For properly chosen dot and material 



parameters, carriers become confined in all three spatial directions within the low- 
bandgap islands on a typical length scale of tens of nanometers. This three-dimensional 
confinement results in atomic-like carrier states with discrete energy levels. When 
the semiconductor is optically excited, an electron is promoted from a valence to 
a conduction band. In the usual language of semiconductor physics this process is 



described as the creation of an electron- hole pair or exciton [5J; [55]: the electron 
describes the excitation in the conduction band, and the hole accounts for the properties 
of the missing electron in the valence band. Electron and hole are conveniently 
considered as independent particles with different effective masses, which mutually 
interact through the attractive Coulomb interaction. 

In the so-called strong confinement regime the wavefunction extension is comparable 
or smaller than the excitonic Bohr radius of the bulk material, and the electron-hole 
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states acquire a dominant single-particle character with only slight renormalizations 
due to Coulomb interactions. This situation approximately corresponds to that of most 
types of self-assembled quantum dots [H; 0; 0] where carriers are confined in a region 
of typical size 10 x 10 x 5 nm 3 . To the lowest order of approximation, the groundstate 
of the interacting electron-hole system is simply given by the product of electron and 
hole single-particle states of lowest energy <p e (r) and (^(r), respectively. Within this 
work we shall consider the situation where an empty quantum dot is optically excited. 
By tuning the frequency to the groundstate-exciton frequency and choosing the proper 
light polarization, fig. [TJ one can selectively excite one of the two spin-degenerate exciton 
states of lowest energy. The population of excited exciton states can be completely 
suppressed owing to their large energy separations of several tens of meV. Alternatively, 
one can also excite within a two-photon process a biexciton consisting of two electron- 
hole pairs with opposite spin orientations. In comparison to the exciton energy uq, 
the biexciton energy 2 ujq — A is reduced by the biexciton binding A ~ 2 meV due to 
Coulomb correlation effects, which allows to individually address excitons and biexcitons 
through proper frequency selection. 

In this work we shall investigate optical experiments with short laser pulses, where 
the semiconductor vacuum is coupled to only the exciton or biexciton groundstate and 
population of other excitonic states is suppressed because of the large energy splittings. 
Introducing the usual notation of for the vacuum and 1 for the exciton or biexciton 
state^the system's dynamic in presence of light coupling is governed by the Hamiltonian 



53; 56 



H = E |1)(1| + (e-^(t) |1)(0| + e^Q*{t) |0)(1|) , (2.5) 

with Eq the exciton or biexciton groundstate energy, uj the central frequency of the laser 
pulse, and Q(t) the Rabi frequency associated to the laser pulse envelope. The latter 
quantity determines the effective light-exciton coupling, and depends on the electric field 
strength of the laser pulse and the dipole moment of the excitonic transition. In deriving 



eq. (12.51) we have made use of the usual rotating- wave approximation [561 ] . The coupling 
of the exciton or biexciton groundstate to phonons is governed by the matrix elements 
(l\H cp \l). Phonon couplings to other excitonic states can be safely neglected because of 
the large energy separation from the semiconductor groundstate (gap energy E g ~ eV) 
and from excited excitonic states (several tens of meV), where typical acoustic-phonon 
energies are of the order of only a few meV. The total system Hamiltonian thus reads 

H = H Q + ^u q a\a q + ^g q (a q + a) q ) |1)(1|. (2.6) 

Q Q 

Here u) q = c^q is the acoustic phonon energy with q the semiconductor sound velocity, 
and g q the exciton-phonon matrix element assumed to be real-valued. For the product- 
type exciton wavefunction e (r)0h(r) considered above, the matrix element reads 



28; 31 



9 9 =[^) [ D e I dre~ iq ' r \<p e (r)\ 2 - D h \ dr e~ iqr \(j) h (r)\ z ) . (2.7) 
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The integration terms on the right-hand side are precisely the form factors for electrons 
and holes, respectively. For a biexciton consisting of two electron-hole pairs with 
opposite spin orientations, the matrix elements (12.71) have to be multiplied by an 
additional factor of two. 



3. Analytic results 

The system described by the Hamiltonian (12.61) is known as the independent Boson 
model 28]. In contrast to the phonon couplings (12.41) . the independent Boson 
hamiltonian does not induce transitions between different states. Yet it leads to 
decoherence. This can be easily seen by writing eq. H 2 . 6 1) in the interaction picture 
according to the hamiltonian for free phonons, viz. 

V(t) = (e-^ f a q + e^a\) |1)(1| . (3.1) 

Q 

Through the phonon coupling the quantum dot state becomes entangled with the 
phonons, where each phonon mode evolves with a different frequency u q . If we trace 
out the phonon degrees of freedom, as discussed in more detail below, the different 
exponentials e * interfere destructively, which leads to decoherence. Because this 
decoherence is not accompanied by relaxation, the process has been given the name 
pure dephasing. We shall now study things more thoroughly. Suppose that an ultrashort 
laser pulse brings the system at time zero from the groundstate into the superposition 
state 

if) = cos#|0) + sin0|l) , (3.2) 

with 9 a small mixing angle. We shall study next how ip evolves in presence of phonon 
couplings. In doing so we make the reasonable assumption that the phonons are initially 
in thermal equilibrium and decoupled from the quantum dot, which will allow us to 



solve the problem analytically [28|; 141]. To understand the essentials of this scheme, 
it suffices to consider a single phonon mode. The response to the phonon continuum 
results, as suggested by the term "independent-Boson model", by simply summing up 
the contributions of all phonon modes. 



3.1. Single phonon mode 
Let 

h = e |l)(l| + u a^a + g (a + a ] ) |1)(1| (3.3) 

be the Hamiltonian of a single phonon mode interacting with a two-level system. We 
next introduce the phonon displacement operator 56|; |57|; [58 

J D(0 = exp(ea t -r«), 



(3.4) 
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which has the properties Dl(£)aD(£) = a + £ and £>t(f) t£)(£) = a t + f*. With the 
operator s = a{a^ — a)|l)(l| one then finds e s = |0)(0| + where a = g/u. 

Thus, 

h = e~ s (e(,|l)(l| +uo t a) e s , (3.5) 

with Cq = tQ—g 2 /uj the renormalized energy of the two-level system. In other words, the 
independent Boson hamiltonian (13.31) results from the hamiltonian ho = e' \l)(l \ +ua f a 
by shifting through D(a) the phonons when the system is in the upper state, accounting 
for the fact that the exciton provides a perturbation to the lattice which becomes 
distorted. This is reminiscent of molecular physics where electronic excitations are 
accompanied by a change of the binding properties and consequently the molecular 
structure, though in our case the coupling is much weaker and to a continuum of phonon 
modes rather than to a few vibronic states. Consider next the time evolution operator 
in the interaction picture 

U (t, t') = e ihat e -ih(t-t') e -ih t' = e ih t e -s e -ih t e ih t' e s e -ih t' ^ (3 ^ 

where we have used eq. (13.51) to arrive at the final expression. The first three terms of 
the last expression can be simplified to 

e ih t e -s e -ihot = | > <0| + D(-ae^*)|l)(l| , (3.7) 

and a similar expression follows for the last three terms. Using the property £)(£)£)(£') = 
-0(£ + £')e* 3m ^' we then arrive at our final expression 

u (t,0) = |0)(0| + e- ia2sin ^ J D(a[l-e^'])|l)(l|, (3.8) 

where we have set for simplicity if — 0. 

3.2. Phonon continuum 

It turns out that a completely similar procedure can be applied in case of a phonon 
continuum. This is because each phonon mode couples independently to the exciton or 
biexciton, and the response of the system can be simply computed by summing over 
different modes. With s = J2 q a q( a q ~ a 9 )|l)(l| an d ct q = g q /u) q one then finds, in 
analogy to eq. ( 13. 5ft . 

/J = e-(^|l)(l| + ^^ata q )e s , (3.9) 

Q 

with Eq = Eq — ^2 q gq/w q the renormalized exciton or biexciton energy. For the time 
evolution operator one obtains after some straightforward calculation 

U(t,0) = \0)(0\+Y[(e- i < s ™"« t D q (a q [l-e i "« t })) |1)(1|. (3.10) 
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Table 1. Dot and GaAs-type 


material parameters used 


in our calculations. For the 


electron and hole wavefunctions e (r) and 4>h{r) we assume Gaussians with full width 


of half maximum (FWHM) Li 


and L z , respectively. 




Property 


Symbol 


Value 


Mass density 


P 


5.37 g/cm~ 3 


Sound velocity 


C£ 


5110 m/s 


Deformation potential for electrons 


D e 


-14.6 eV 


Deformation potential for holes 


D h 


-4.8 eV 


in-plane confinement (FWHM) 


U 


8 nm 


confinement in z-direction (FWHM) 


L z 


3.5 nm 


biexciton binding 


A 


2.75 meV 



3.3. Dephasing 

We next consider the situation where initially the quantum dot is in the groundstate and 
the phonons in thermal equilibrium, and the quantum dot is brought into a superposition 
state (13 .2p at time zero. Correspondingly, the phonons will react at later time to the 
perturbation of the excitonic system. We shall only be interested in the quantum dot 
system and hence trace over the phonon degrees of freedom. To account for this open- 
system dynamics we employ a density-operator framework. Let p = ® P P h 
be the initial density operator, with ip the initial superposition state fl3.2l> and p p h the 



operator for phonons in thermal equilibrium 281 1 . Owing to the electron-phonon coupling 



(12.61) . which preserves the exciton number, the lower- and upper-state populations 
Poo(t) — c °s 2 an d Pn(t) = sin 2 9 do not change with time. In contrast, the polarization 

p w (t) = (l\tr ph U(t,0)p o U(0,t)\0) (3.11) 

does. To compute this polarization function, we use the time evolution operator (13.101) 
and utilize that (0|C/(0, t)|0) = 1. Then, 

pio(i) = sin^ cos 9 tr ph JJ e~< A ^D q {a q [\ - e^*]) . (3.12) 



It can be shown [28J; |59|; [60] that in thermal equilibrium the expectation operator 
of the displacement operator is (_D(£)) = e _ ^' 2< - n+ 2\ with n the usual Bose-Einstein 
distribution function for a given phonon mode. Putting together all results and 



performing some simple manipulations we arrive at our final result 53 



/,\ • n r\ ( S3. \ (i s'mu>nt+(l— cos u)„t) coth ^2.^) /„ -, \ 

p 10 (t) = sm^cos^e *-">\"<iJ \ q y q> 2 J, (3.13) 

with (3 the inverse temperature. Equation (13. 13[) is exact and holds for arbitrary coupling 
constants g q . In evaluating the sum in the exponent one replaces — > J d 3 q/(2ir) 3 to 
account for the continuum of phonon states. 

Figure [2](a) shows the correlation function p 10 (t)/(sin6'cose) for the realistic dot 
and material parameters listed in table [TJ and for zero temperature. The correlation 
function drops on a picosecond timescale to a value of approximately 0.9 and 0.6 for 
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Figure 2. (a) Correlation function (|3. 13|) for exciton X and biexciton XX, and for 
the dot and material parameters listed in table [TJ (b) Density plot of the phonon 
displacement (|2.3|) as a function of time and position, and for the situation where 
the system is brought at time zero to the superposition state ()3.2|) . The dashed line 
indicates the extension of the exciton wavcfunction, and the arrow points to the emitted 
phonon wave-packet. Panels (c-e) report snapshots of the phonon displacement for 
selected times of 0.1. 1, and 2 ps. 



excitons and biexcitons, respectively. This coherence loss is due to phonon couplings, 
and is stronger for biexcitons than for excitons because of the larger coupling constant. 
The picosecond timescale is determined by the maximum of the electron-phonon 
coupling constant (12. 7p . which, in turn, is governed by the electron and hole form 
factors. Thus, for small dots the correlation drops faster than for large dots because of 
the larger portion of phonon modes to which carriers in small dots can couple. Phonon- 
assisted dephasing was first observed by Borri and coworkers 24j; |25|; |26j in four-wave 
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mixing experiments, and was confirmed in a l arg e number of independent studies. 
It constitutes a major obstacle for optical Il5t Il6| and electrical [48| quantum- 
information applications as it spoils the ideal quantum state evolution. 



3. 4- Where does the quantum information go? 



As the time evolution described by the operator. ( 13.101) is completely unitary, one might 
wonder how the quantum coherence becomes reduced. In fact, at zero temperature the 
composite dot-phonon system is described by a single wavefunction. As time goes on 
the dot degrees of freedom, i.e., exciton or biexciton, become entangled with the phonon 
degrees of freedom. It is thus the tracing procedure in (13.111) that provides decoherence 
of the quantum dot system. Here our ignorance about the environment, i.e., phonon 

Where does the quantum 



61 



degrees of freedom results in a partial coherence loss 
information go? To answer this question we plot in panels (c-d) of fig. [2] snapshots 
of the lattice displacements (12.31) subsequent to the preparation of the superposition 
state at time zero. One observes that it takes some time for the lattice to react 
to the excitonic perturbation. This is shown even more clearly in panel (b) which 
reports the time evolution of the lattice distortion along a given direction x. Due to 
the phonon inertia the lattice distortion overshoots, instead of smoothly approaching 
the new eq uilibrium position, and a phonon wave packet is emitted from the dots (see 



arrow) [35j; |48|; [62| which imprints the quantum information about the superposition 



state into the environment and thus reduces its quantum properties: the system suffers 
decoherence. It is thus primarily the abrupt quantum-state preparation that causes 
excitonic decoherence. We will return to this point later when discussing how such 
decoherence can be efficiently suppressed by means of quantum control. 

4. Wavefunction approach 



Equation (I3.13P is an exact result that was derived for a delta-like optical excitation. 
No comparable simple expression can be found for optical excitations with laser pulses 
of finite temporal width. To account for such situations, one can either approximate 
the pulse by a comb of delta pulses 32|; [33( or employ a density-matrix description 



36 



47; 53 



In this work we shall use the latter approach. For simplicity we shall 
restrict ourselves to zero temperature where the composite dot-phonon system can be 
described by a single wavefunction. The treatment of finite temperatures within a 



density- matrix approach is similar but more tedious [36|; |47J; |53|. Quite generally, our 



approach is expected to be appropriate also for finite temperatures that are sufficiently 
small in comparison to those phonon energies cu gmax where the electron-phonon coupling 
fiWx i s maximal. For the dot parameters listed in tabled] we find uj qmaK ~ 10 K. 

To simplify our analysis, we shall introduce some approximations that are motivated 
by the analytic expression (13.131) . We first consider an isotropic quantum dot model 
with an effective phonon coupling g q = (q q )n which is obtained by averaging g q over 
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all angles. One can easily check that the correlation function (I3.13P is not altered by 
this averaging procedure. We next approximate the phonon continuum by a finite set of 
effective phonon modes. Let p denote an index that labels the modes with wavevectors 
q p = pAq. Here Aq is the wavevector spacing. To account for the three-dimensional 
phase space of phonon modes, we have to modify the effective phonon coupling according 
to g p = (q p Aq/ (27r) 3 )5 g Qp to arrive at our effective Hamiltonian 

H = H Q + Y^ UJ P h l% + H~9p{' 1 P + ~ a l) |1><1I , (4-1) 
p p 

with H the quantum-dot hamiltonian (12.51) and ap the creation operator for the effective 
phonon mode p. We next specify our computational Hilbert space. Starting from the 
quantum dot states |0) and |1) we create one-phonon excitations through al\0) and a p \l), 
two-phonon excitations through a p a p ,\ 0) and a p a p/ 11), and so forth. The truncation of 
this Hilbert space basis is performed on the basis of eq. (13.1 3p . which shows that the 
importance of phonon modes to the correlation function is governed by a q = g q /oj q . We 
correspondingly keep in our Hilbert space only one-phonon excitations with a p > e a max , 
two-phonon excitations with a p a p > > ea max , and so forth. Here a p = g p /uj p , a max is the 
maximum of the a coefficients, and e is a small cutoff parameter. In our calculations 
we typically keep hundred phonon modes and use a cutoff parameter of e = 10~ 2 . 



4-1. Comparison with analytic results 

The dotted lines in fig. [2] (almost indistinguishable from the solid and dashed lines) show 
results of our wavefunction simulations for the situation where the dot is initially in the 
superposition state (13. 2p and the phonons unexcited. One clearly observes that the 
results of this approximate wavefunction approach is in perfect agreement with those of 
the analytic expression (13.131) . thus proving the accuracy of our scheme. We emphasize 
that the cutoff based on the a values is crucial to obtain such nice agreement, whereas 
a mere restriction to one- or two-phonon processes would not work equally well. 



4-2. Rabi oscillations 



We now turn to the situation where a quantum dot initially in the groundstate is excited 
by a short laser pulse. Throughout we use Gaussian laser pulses with a full width of 
half maximum of 2.33 ps, approximately corresponding to the experiment of Zrenner 
et al. 



18 



21 



63 



and assume that the laser frequency is tuned to the renormalized 
exciton or biexciton transition energy E' , respectively. Figure [3] shows results of our 
simulations. We shall find it convenient to analyze our results in terms of the usual 
Bloch vector u 53j; |54J |: the z component accounts for the population inversion which is 
— 1 when the quantum dot is in the groundstate and 1 when it is in the excited state; the 
x and y components account for the real and imaginary parts of the quantum coherence 
pio, respectively. The length of the Bloch vector is one for a pure state, and smaller 
than one when the system has suffered decoherence. 
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Figure 3. Rabi oscillations in presence of a Gaussian laser pulse for: (a) exciton 
and pulse area 9 — 2n, (b) biexciton and 9 = 2?t, and (c,d) exciton and biexciton for 
9 = 4ir. The solid line shows the z-component of the Bloch vector and the dashed line 
its modulus. The insets report the Bloch vector trajectories. 

Panel (a) in fig. [3] shows the Rabi oscillation of an exciton for a laser pulse with area 
9 = £l(t)dt = 2n. The inset reports the trajectory of the Bloch vector during the 
laser excitation. The solid line shows the time evolution of the ^-component of the Bloch 
vector, which passes from the groundstate through the excited state (exciton) at time 
zero back towards the groundstate. However, because of phonon couplings the systems 
loses part of its coherence to the phonons, and as consequence does not fully return 
to the groundstate, as it would in absence of phonon couplings. Instead it remains in 
an excited state: the system has suffered decoherence. This is also apparent from the 
dashed line that reports the length of the Bloch vector, which decreases from a value of 
one to approximately 0.9, thus indicating the transition from a pure to a mixed state. 
Similar behavior is observed in panel (c) for a pulse area 9 = 4ir, where the Bloch vector 
is rotated twice. Finally, panels (b) and (d) show simulations for Rabi floppings of 
biexcitons, which couple more strongly to the phonons. Correspondingly, the coherence 
losses are more pronounced. In particular for the 4ir rotation shown in panel (d) the 
system ends up in a mixed state with equal probabilities for the the upper and lower 
state populations. We note that these excited exciton or biexciton states would decay 
on a much longer time scale of nanoseconds due to radiative recombination. 

In fig. H] we plot the upper-state population p n as a function of pulse area. For 
the excitons (solid line) one clearly observes Rabi-type oscillations which are damped 
for larger pulse areas. Note that for all pulse areas 9 we have used the same pulse 
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Pulse area 9 (n) 

Figure 4. Upper-state populations after completion of Rabi oscillations as a function 
of pulse area. 



widths. Thus, for large ^-values the system oscillates very fast between ground and 
excited state, and becomes dynamically decoupled from the phonons (see also discussion 
further below). This is reflected by the fact that the oscillations in fig. H] are less damped 
for large 9 values. For the biexcitons we use, contrary to the excitation scenario used 
above, i.e. a Gaussian pulse with given area, a situation that is related more closely to 
experiment. In fact, in a two-photon experiment the effective Rabi frequency is given 
by f2 c fj = fi 2 /(2A) 6J; |65(, with A the biexciton binding energy (see tabled]). The 
corresponding Rabi oscillations are shown by the dashed line in fig. HI One observes a 
different period of the oscillations and a substantially stronger damping in comparison 
to excitons. We finally emphasize that our results are in almost perfect quantitative 



agreement with the experimental results of refs. 2l|; |65 



5. Quantum control 

Above we showed for the independent-boson model that in presence of an exciting laser 
pulse Rabi flopping occurs but is damped because of phonon decoherence. However, 
contrary to other decoherence channels in solids where the system's wavefunction 
acquires an uncontrollable phase through environment coupling, in the independent 
Boson model the loss of phase coherence is due to the coupling of the electron-hole 
state to an ensemble of harmonic oscillators which all evolve with a coherent time 
evolution but different phase. This results in destructive interference and dephasing, 
and thus spoils the direct applicability of coherent carrier control. On the other hand, 
the coherent nature of the state-vector evolution suggests that more refined control 
strategies might allow to suppress dephasing losses. 



Quantum control of polaron states in semiconductor quantum dots 



14 



5.1. Quantum bang-bang control 



A first step in this direction was undertaken by Viola and Lloyd 4]J who considered the 
situation where a two level system is initially prepared in the superposition state ( 13.21) 
and at later times a sequence of 7r-pulses is applied. The operator accounting for tt- 
rotations is U n = |0)(1| + |1)(0|. Similar to the derivation of the time evolution operator 
I3.10p . one can show that the time evolution for a sequence of two such tt- rotations is 



41] 



U(t + 2St, t + St)U n U(t + St, t)U n 



|0> <0| + |1)<1|JJD, (-e^'a, (e^-l) 2 ) ~ 1 + 0{5t 2 ) 



(5.1) 



where pure phase factors have been suppressed. The important observation of Viola 
and Lloyd concerns the fact that for sufficiently small time intervals St, i.e., much 
smaller than the time on which quantum coherence is lost, the system does not suffer 
decoherence. Rather its dynamics is governed by fast Rabi flopping such that it becomes 
decoupled from the environment. Consequently, this control strategy has been given the 
name "quantum bang-bang control" . 



5.2. Optimal quantum control 

One might wonder whether there exist more simple control strategies that allow to fight 
decoherence. A general framework for the determination of such control is provided 



44 



45 



46 



53J. Here, the objective of the control is 



by optimal control theory (OCT) 
quantified through a cost function, which is then minimized subject to the condition that 
the time dynamics of the coupled dot-phonon system is governed by the Schrodinger 
equation. Optimal control theory is a mathematical device that has found widespread 
applications for, e.g., molecules 



45 



66 



67|, or semiconductors |47l ; 



48 



atoms 

In the following we shall analyze the situation where a 7r-pulse brings the quantum 
dot from the ground to the exciton or biexciton state. Our control objective is to achieve 
a maximum upper-state population after completion of the pulse. In the upper state the 
density operator is of the form |l)(l|®p p h, i-e., the quantum dot is in a pure state. Thus, 
the transition from the ground to the excited state can only be achieved if decoherence 
losses are fully suppressed. Let T denote the pulse duration and ip(T) the terminal 
wavef unction. The cost function 

1 



M1>) = MT)\F \iP(T)) 



(5.2) 



rates the success of a given control fl(t). Here P = |0)(0| <8> l p h is the projector on 
the quantum dot groundstate. Jt is minimal when the system ends in the upper state, 
irrespective of the final phonon state. To make our optimal-control problem well posed 
we add an additional term to the cost function, viz., 



7 



n{t) 



dt 



(5.3) 
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that favours control fields Q(t) with a smooth time variation. 7 is a weighting parameter 
that determines the importance of the two different control strategies of upper-state 
population and smooth control fields. We shall use small 7 values throughout, such that 
the cost function J(if>, Q) is dominated by the first term. The control problem under 
consideration thus becomes the minimization of the cost function J(ip, fl) subject to the 
condition that ip(t) fulfills the Schrodinger equation. Within the field of optimal control 
theory one uses Lagrange multipliers to turn this constrained minimization problem into 
an unconstrained one. For this purpose we define the Lagrange function 

L(i>,p, n) = J(V, fi) + (p(t), ij>(t) -(H + H op (t)) m) , (5-4) 

with H the independent Boson hamiltonian (12.61) . H op the light-matter coupling ( 12. 51) . 
(u,v) = Jq dt(u\v), and p(t) the Lagrange multiplier. We next utilize the fact that 
the Lagrange function has a saddle point at the minimum of J(ip,Q), i.e. all three 
derivatives 5L/5ijj, 5L/5p, and 5L/5Q must be zero. Performing usual functional 
derivatives in eq. JS7 



optimality system [46 



1) we obtain after some variational calculation the following 



53 



iif> = {H + H op (t))iP (5.5a) 
ip = {H + H op {t))p (5.55) 

7 fi = -Ke(^|(|0)(l|®U ph )b) ; (5-5c) 
which has to be solved together with the initial and terminal conditions 

V(0) = ^0 (5.6a) 

ip(T) = F ip(T) (5.6b) 

fi(0) = fi(r) = 1 . (5.6c) 

The right-hand side of eq. (15.6fcj) follows from the functional derivative 5 J/ Sip. Notice 
that while the state equation (15.5oj) with initial condition ^(0) = ipo evolves forward in 
time, the adjoint equation (I5.5fc|) with terminal condition (I5.6fc|) is marching backwards. 
The control equation (I5.5cl) determines the optimal control. In most cases of interest 
one is not able to directly guess Q(t) such that eqs. (l5.5aL c) are simultaneously fulfilled, 



and one has to employ an iterative scheme. In this work we follow ref. [46j and 
formulate a numerical algorithm that solves the optimality system (15.5 d -c) for a given 
initial configuration tp Q . To solve this problem, we apply a gradient-type minimization 
algorithm, which, starting from an initial guess for Q(t), determines a search direction for 
an improved control. We first solve eq. (I5.5al) with initial condition ^(0) = ipQ forwards 
in time. Once the wavefunction if>(T) at time T is computed, the final condition for p(T) 
can be calculated from eq. ( 15.661) and the adjoint equation of motion ( 15.561) is solved 
backwards in time. The gradient of L with respect to Q becomes 

_ = _ 7 a - Ke^KVl (|0)(1| ® l ph ) \p) , (5.7) 

which gives the search direction for an improved control that minimizes J(i/j,Q). In 
the following we employ for the minimum search the usual nonlinear conjugate gradient 



method 68 
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Figure 5. Results of our calculations with a Gaussian it (gray lines) and OCT (dark 
lines) laser pulses. Panel (a) shows f2(t) and panel (c) the time evolution of the z- 
component (solid line) and modulus (dashed line) of the Bloch vector for excitons. 
The insets report the trajectories of the Bloch vector for the OCT (left) and Gaussian 
(right) pulse. For the Gaussian 7r-pulse Rabi flopping occurs but is damped due to 
electron-phonon interactions. For the optimal control decoherence losses are completely 
suppressed, and the the system ends up in the desired upper state at time T. (b,d) 
Same as (a,c) but for biexcitons. 



Results of our optimal-control calculations are shown in fig. [51 Panel (a) reports 
the Gaussian (gray line) and optimized (dark lines) control fields Q(t) for excitons. The 
corresponding Bloch vector trajectories are depicted in panel (c). We first observe 
that for the Gaussian control field Rabi flopping occurs, but the system ends up 
in a mixed state owing to phonon-assisted decoherence. In contrast, the optimized 
^oct(^) channels the system from the ground to the excited state without suffering 
any decoherence losses. Thus, optimal control theory allows to design control strategies 
that can drastically outperform more simple schemes. In ref. |4Sj we showed for a related 
system that the optimized quantum field strongly suppresses the emission of a phonon 
wave packet, and allows the lattice to react smoothly to the time varying excitonic 
configuration. In contrast, for the delta-like control shown in fig. [2] as well as the 
Gaussian pulse the superposition properties of the quantum dot state are imprinted 
into the lattice environment by emitting a phonon wave packet, and the system suffers 
decoherence. Similar conclusions apply for the biexciton control shown in panels (b) 
and (d). 

Our optimal quantum control strategy differs appreciably from other control 
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strategies. The inherent coupling of electrons_to phonons excludes quantum state 

or other quantum-optical control 



42 



manipulations in decoherence free subspaces 
techniques, such as, e.g., stimulated Raman adiabatic passage [69(, where quantum state 
transfer is achieved through states fully decoupled from the environment. Furthermore, 
the time dynamics of the phonon degrees of freedom disables spin-echo techniques to 
restore pure quantum states by means of effective time reversal through n pulses. Finally, 
optimal control only requires smooth laser pulses on the timescale of picoseconds, 
rather than femtosecond pulses needed for quantum bang-bang control [41] . where the 
system has to become dynamically decoupled from the environment. We note that an 
alternative to laser-pulse shaping might be to use laser pulses with fixed, e.g., Gaussian 
shape, and to perform the quantum control by means of voltage pulses applied to 
an external gate. Here, the energy detuning between the central laser frequency and 



the two-level system can be modified .65]. Corresponding optimal control calculations 
indeed demonstrate the applicability of optimized voltage pulses for the purpose of 
efficient exciton flopping in presence of phonon-assisted dephasing. 



6. Conclusions and summary 

In conclusion, we have studied the phonon- induced dephasing dynamics in optically 
excited semiconductor quantum dots within the frameworks of the independent Boson 
model and optimal control. We have shown that appropriate tailoring of laser pulses 
allows to control the dot states without suffering significant dephasing losses. The 
requirements for such laser-pulse shaping are well within the possibilities of presentday 
technology. We attribute our finding to the fact that in the process of decoherence it 
takes some time for the system to become entangled with its environment. If during 
this entanglement buildup the system is acted upon by an appropriately designed 
control, it becomes possible to channel back quantum coherence from the environment 
to the system. We therefore believe that our findings are relevant for a much broader 
class of solid state systems where quantum information is encoded in long-lived quasi- 
groundstates with small energy separations, such as electron or nuclear spins, resulting 
in slow scattering processes that can be manipulated by means of quantum control. 
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